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1 Introduction 

Since it was first proposed by Brandenberger and Vafa [H [2 [31 H] , the string gas cosmology scenario has 
generated a significant amount of interest. One of its most appeahng characteristics is that it provides 
a mechanism for dynamically generating a four dimensional space-time. This argument is based on the 
assumption that strings interact mainly by intersecting each other. If that is the case, the probability of 
intersection in space-time of two worldsheets has non-zero measure only if the dimension is equal or less than 
4. This is a classical argument and it is not obvious that it will remain true if quantum effects are taken 
into account. There have been several attempts at trying to formulate and prove the Brandenberger- Vafa 
mechanism with mixed results O [6l [7l [HI |9j . A recent work [10], for example, succeeds in decompactifying 
3 large spacial dimensions for a gas of diluted strings. 

Apart from the Brandenberger- Vafa mechanism, there have been other attempts to produce a mechanism 
for realizing a four dimensional space-time |111 112| . One of these scenarios consists in the inclusion of a 
two-form field. This field is already present in the supergravity action, hence, it is natural to consider 
its appearance in the equations of motion. Cosmologies with a two-form field had been studied in the 
past and several solutions are known |13j . In the context of string gas cosmology, a two- form field flux 
was introduced for dilaton-gravity in [HI [15]. In these solutions the 2-form field flux is restricted to a 
4-dimensional submanifold of space-time. The two-form field flux introduces a repulsive potential in the 
equations of motion for the scale factor of the spacial dimensions. As a result the expansion of three spacial 
dimensions is enhanced and the corresponding scale factors become large. 

In [16 , Greene et al. introduced a higher derivative dilaton gravity model. This model replaces the 
Newtonian-like kinetic terms in the dilaton gravity action by their relativistic counterparts. By doing this, 
one obtains a model with some nice features: derivatives with respect to the cosmic time become bounded, 
singularities at finite time are avoided, bounces on the scale factor are produced and loitering phases that 
solve the horizon problem are realized. Since this model respects T-duality, it has been applied to investigate 
stringy cosmology such as the Brandenberger- Vafa scenario. The results has shown that this model could 
lead to a 4 dimensional space-time only with a fine tuning. In general 3 large spacial dimensions are not 
preferred and any number of dimensions could become large. 

Because of the many appealing features of this model, it is interesting to investigate its behavior in the 
presence of an antisymmetric tensor field. The effects of such kind of field have been studied in several works 
[T51 FT7| in the context of dilaton gravity. There have been some studies also for the string gas cosmology 
[m [TS] case. For string gas cosmology, it is specially interesting to introduce a 2-form gauge field. Since the 
initial configuration is supposed to be compactified on a d-dimensional torus, there exist non-trivial effects 
even for a constant gauge field as strings can wrap the compact dimensions. 

In section 2, we briefly recall the low-energy string effective action. In section 3, we introduce a T-duality 
invariant effective action with 2-form flcld flux. We obtain analytic solutions of the effective action in several 
limiting cases. We find they can explain qualitative behaviors of the numerical solutions. One of the main 
features of the string cosmology model is the introduction of a Hagedorn phase for the early universe. It 
arguably removes the initial singularity of the universe. We find that a homogeneous 2-form field flux triggers 
an expansion of a universe away from the oscillating phase around the self-dual radius. An interesting stringy 



effect, as explained in section |4j is that an electric like two- form field modifies the effective string tension 
and the Hagedorn temperature [HI [19]. In such a case, the energy of the winding modes can vanish for the 
directions parallel to the electric field. We find that these spacial directions can expand even if the winding 
modes are present. We conclude in section 5. 

2 Low-energy string effective action 

In this section we recall the low energy effective action for string theory. 
The low energy effective action in the string frame is given by |20| 

^ = ^ y d^\V^e-^ (r + (V0)2 -V- ^hA + I d^\^^C,natter (2.1) 

where k^q = SnGio and H^ — Hfj^^^H'^^^ with if^i^A = d[f^B^x]. The sign convention is all + according to 
the classification in Misner, Thorne and Wheeler [3T]. The variation of this action gives the equations of 
motion 

- Ig^i^^r + igf.g"' - 9yWxy.<t> (2.2) 

W^,{e-*H^"'^) = (2.3) 

2n<j) + R-[\/(j)f ~V - ^H'^ = Q (2.4) 

where T^i, is the energy-momentum tensor derived from the matter Lagrangian. We assume the space-time 
metric is of the following type in the string frame 

ds^ = -df + a^dxidx' + b^dxidx' , i = {1, 2, 3} , / = {4, ..., 9} (2.5) 

with 

a = e^(*), b = e''^'^ (2.6) 

By a conformal rescaling 

5^1^ = e^^Ofiiy (2.7) 

we obtain the effective action in the Einstein frame 



+ / d^° Xy/^ e^^^Cmatter (2.8) 

The field equations are 

Rf^u - \gt.uR = k?o(T'm. + ''"^T^u + '^^^T^. + ^^^7;.) (2.9) 



V^(e-^iJM>^-^)^0 (2.10) 

U4>-Vei + -e^'^H^ = Q (2.11) 
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The homogeneous metric is given by 

dS^ = -dP + d^dx.dx' + b^dxidx' , i = {1, 2, 3} , / = {4, ..., 9} (2.12) 

with 

a = e"(*), 6=e^(*) (2.13) 

More detailed relations between the string frame and the Einstein frame are explained in the appendix A. 
We need a nontrivial solution for the field strength Hf^^\ to investigate its effects on the cosmology. The 



equation of motion for the two-form field (2.10) can be solved using the Freund-Rubin ansatz 



H^""" = e'f'e>"'"^Vph , with e^"^*^ = -^S\'^d'[6^5^-^ (2.14) 

We assume here that the two- form field flux exists only in three spacial directions. This assumption is 
consistent with the symmetry of the postulated space-time metric. Because V^e^j/Q.^ = 0, the equation of 
motion is automatically satisfied. It only remains to satisfy the closure condition |13) 

V[^i?^,„] = (2.15) 

The equation of motion is then 

3 

^ 5^''a^(a3r6e^V) = (2.16) 

fi, ly—Q 

In the case of a homogeneous field h = h{t), we obtain 

h+{3X-6i> + ^)h = (2.17) 

with Ho a positive constant and ' = ^ • This equation is solved by 

h=^^e-* (2.18) 

Thus, we have for these H^^\ 

H^XnH''^^ = if M,^- {0,4, 5,. ..,9} (2.19) 

H.xnH^^^ = '^Si ioT I, J ^ {1,2,3} (2.20) 

In particular we have H^ = 6H^a~^. 

As in [I] [51 [TH [TC] we consider a very simple setup with 3 types of matter: isotropic winding modes (with 
all winding numbers Wi = W, i = {1, 2,3}) with energies 

Ew = QWe^ (2.21) 

isotropic momentum modes (with all momenta Ki = K, i = {1, 2,3}) with energies 

Ek = GKe"^ (2.22) 



and string oscillator modes that are modeled as pressureless dust with energy -Edust Q The total energy is 
the sum 

E^Ew+Ek + Ednst + V (2.23) 

with V = V{ip) the potential for the dilaton. In an adiabatic system the pressures are 

BE _dF _ dL„, _ dV 

dip dip dip dp 

'^ 

2Ke-^ - 2We^ (2.25) 



P. - ^ = 1^ = -^-^ (2-24) 

1 dF 1 ai,„ idE _, _^ ^^^^ ^ 



3dX 3 dX 3dX 

The energy of the string gas is defined as Eg = Ew + Ek + £^dust • 

In order to model the behavior of the gas, we consider the following phases as in |16j : 

• Hagedorn phase: thermal equilibrium at temperature Th — l/(-\/87r). The free energy of the gas 
vanishes (Pa = 0) and Eg is conserved 

• Radiation phase: thermal equilibrium at T < Th with the universe dominated by massless string 
modes. In d + 1 dimensional space-time, the internal energy is 



E.,.,y,T'*^, ,,,:,8_^g|i_lL„_,-.) 



(2.27) 



with Vd — (27r)'^e''l'^l: the T-duality invariant volume. 



F^Es~TS = -\cdVdT''+^ (2.28) 

a 

Px = sign{X)Es/d (2.29) 



A > : (VF) = , (K) = ^P^e^ (radiation phase) (2.30) 

A<0: (W) = --Pxe^^ , {K) = (winding mode dominated phase) (2.31) 

Note that the radiation and winding mode dominated phases are T-dual to each other. 

• Frozen phase: in this phase the interactions between strings are turned off. The momentum and 
winding numbers are conserved, so K and W are frozen at the values they have on Hagedorn exit. 

• Non-equilibrium phase: In order to model the string gas, we consider a phase in which the the temper- 
ature falls below the Hagedorn temperature. Since we also consider the interactions among the strings, 
the expectation value of the momenta and winding number deviate from their equilibrium values such 
that the pressure of the string gas does not vanish. 



^dust also contains the contributions from strings with momenta and windings along 6 extra dimensions. 



3 T-duality invariant action with two-form field 

In order to analyze the effect of the two-form field, we work in the string frame. This allow us to choose a 



solution where the scale factor v, defined in (2.6), becomes constant and the analysis can be restricted to a 



4-dimensional cosmology in the presence of a two-form field. 



In the case of a homogeneous space-time [H1[T51[T7] . the action (2.1) can be reduced to 

S = I dt[47re"'^(dA2 - ^^ - Uo{X)) + L^] 



(3.1) 



with Ljn'- the matter lagrangian, ip: related to the original dilaton (j> as (p = 2(p — dX . The potential U(X) 
arises due to the nontrivial two- form field strength -ff^i^A , as shown in [T31 E] 



HUa = 12i?oe 



2„-2dA _ 



= 24[/o(A) 



(3.2) 



The parameter d counts the number of spacial dimensions with the homogeneous scale factor A. Although 
our case corresponds to d = 3 such that the space-time is 4-dimensional, we retain d dependence explicitly 
in the equations of motion in order to keep track of the algebra. 

Now we proceed as (TFj and replace the canonical kinetic terms by their higher derivative extensions. 
This leads to a phenomenological action with bounded velocities (f, X. It thus rules out singularities at any 
finite proper time. With this modification we obtain a higher derivative action for the dilaton and scale 
factor which are coupled to a two-form field strength 



S= dt 



Vi-t^' - V 1 - dx^ - ^i^n + ^ 



(3.3) 



L,n = —F is the negative of the matter free energy (of the string gas) and U{X) is a modified potential as 
explained below. 

String gas cosmology model needs to respect T-duality, a fundamental symmetry in string theory origi- 
nating from the existence of the minimal length (string scale) . It is realized as the symmetry between the 



winding and momentum modes in a toroidal compactification. Since (3.2) is not explicitly invariant under 



T-duality, it is necessary to modify this potential in order to realize the symmetry. Such a modification 
allows us to solve the equations of motion near the self-dual radius numerically. An adequate choice is 



1 



1 



t/(A) ^ -ij2(e2A + e-2A)-d ^ __^2(^^3^2A) 



(3.4) 



as (3.4 1 is not singular at A = and it reduces to (3.2) for large A 



Defining the relativistic factors [16] 



1 



V^^W 



7A 



Vl - dX^ 



the equations of motion obtained from the action (|3.3| are 



7a 






7a" 



7a ^) - A 



f/(A)) 
dU 



dX 



1 

1 

8^ 



d'Xe'^Px 



(3.5) 

(3.6) 
(3.7) 



We also need to impose the Hamiltonian constraint 

lip-l\- U{X) = 



87r2 



-Ee^ 



(3.8) 



where E is the energy contained in matter. Notice that in the positive energy region ^^p — ^\ > U{X). 



The pressure in (3.6|, (3.7) for the dilaton and the scale factor are defined in (2.24) and (2.25). Rendering 



equations (3.6), (3.7) into a more manageable form, we obtain 



"P 



X 



^-%'{^x' + UW-^e^Pp) 



(1 - rfA') 



ipX 



1-.' 



IdU 
dl)X 



1 

8^ 



e'^Px 



(3.9) 
(3.10) 



Before trying to find some solutions to the equations of motion, let us examine the equation ( |3.8[ ) in 
order to get some idea of the expected behavior. If we put U{X) on the right side of the equation, we see 
that A is subjected to the effective potential 

1 



Kff(A) = U{X) 



87r2 



.Ee'f 



(3.11) 



A schematic plot of Voff(A) is presented in figure [T] If we assume, just for the moment, that the dilaton 
has some fixed value, we can observe the dependence of this potential on A. The E dependent term in 



(3.11) grows exponentially as A increases if the winding modes are present. Consequently, this term tends to 



confine the scale factor near the self-dual radius. On the other hand U{X) is a repulsive potential that has 
its maximun value at the self-dual radius. It decreases exponentially as A increases. The term containing 
E, the energy of the string gas, is at the same time modulated by the exponential of the dilaton. Then, as 
(f -> — oo, Vcs{X) flattens for large A. As the confining effect of Vcs{X) diminishes in such a situation, U{X) 
becomes dominant and the scale factor is able to continue growing. It is also possible, depending on the 
initial conditions, for A to undergo oscillations around one of the minima of the potential or the self dual 
radius. In general, as the dilaton is going to weak coupling, these oscillations stop and the scale factor is 
forced to expand by U{X). 

We have explained that the phenomenological action (3.3) with d = 3 is valid for a special class of 



solutions in superstring theory. In this paper, we investigate these solutions in which only the scale factors 
for 3 spacial dimensions are time dependent. The dynamics of this kind of cosmology has been studied before 
in several works [131 [IZl [2U] where solutions for the case d — 3 have been found. 

We should also be careful to point out that physical interpretation may depend on a chosen conformal 
frame. Unless we are able to fix the value of the dilaton cj), we can not consistently conclude that the size 
of a dimension would remain small in the Einstein frame even if it becomes constant in the string frame. 
Nevertheless we may argue that the string frame is theoretically preferred to measure the size of the universe 
as T-duality holds in the string frame. Even with this limitation in mind, we will go ahead to study the 
cosmology in the string frame [5l[10l|14l[T5l|23l|24l[25l|26]. 

3.1 Analytic solutions (d = 3) 

Now we present some analytic solutions that can be obtained by solving the equations of motion. First, we 
assume a simple equation of state P\ = wE, with w a constant and P^p — (no dilaton potential). Using 



X 



Figure 1: Effective potential V^(A) for the scale factor. 



(3.8 1 and the equations of motion we get 
^ = (l-(p2) 
A = (1 - dX^ 



m 



l-{l-^'Y''[{l-dX'f'^-^^{cos\.2X) 



-d 



^^ + 7a ' ( ||(cosh2A)-'^ tanh2A + w{j^ - 7a - -^(cosh2A)-'*) 



2<i 



2M 



(3.12) 
(3.13) 



For the equation of state, we have three specific cases of interest: w = 0, w = \/d and w — —l/d, that 
correspond to pressureless matter, radiation dominated era and winding mode dominated era respectively. 
As the boundary condition for late time asymptotic behavior, we consider 



A^O, (yj— >0, |A|— >oo, (/?— > — oo 



In this limit, the equations of motion (3.12), (3.131 can be approximated as 

w 



A 



^X + i/^g-SrfA 



(ifi^-dX'^ 



Hle-^"^) 



1 



1 



.V' 



-dX 



\2 



ff2 

-"O -2rfA 

2 



3.1.1 Ho = 0, w 7^ case 

We start with the standard string gas cosmology without 2-form field flux. For the case when Hq 
assume the following ansatz 

(f ^ A\ogt + B 
X = Clogt + D 



After substituting them in (3.15) and (3.16), we find 



f = 
X = 



2 

1 + dw'^ 
2w 



1 + dw^ 



\ogt + B 
lost + D 



(3.14) 

(3.15) 
(3.16) 

: 0, we 

(3.17) 
(3.18) 

(3.19) 
(3.20) 



This asymptotic solutions can be seen in figure[2]for d — 3. We have plotted in the same figure the numerical 
solutions for the full equations of motion ( 3.12[ ), ( |3.13 ) with Ho = 0: momentum mode dominated universe 
(green line, w — 1/3), dust dominated universe (blue line, w = 0) and winding mode dominated universe 
(red line, w — —1/3). Of course the green and red lines are T dual to each other. The numerical solutions 
tend to the late time analytic solutions, which are plotted in figure [2] as gray dotted lines. From this solution 
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Figure 2: behavior of A and ip with Hq 
to the late time analytic solutions. 



and w = 0, w 



w 



The gray dotted lines correspond 



we note that the scale factor goes to a constant value very quickly in the absence of any driving pressure. 
This behavior can be seen also in figures [2) [Sl |4) The blue line in every picture represents the case when the 



effect of the pressure and the two-form field vanish, leading to the solution (3.19), (3.20) with w = 0. Notice 
that this solution is valid for arbitrary d. 

For completeness, we mention that there is an additional solution when d = 1, w ^ 0. In this special case, 



assuming ip — —A, the equations of motion (3.15), (3.161 reduce to a differential equation in one variable 
(fi — (^^, (A = — A^). Then we get the solution 



<f = - \og{t + A) + B 
A = \og{t + A) + C 



(3.21) 
(3.22) 



with A, B, C constants. This solution is not physically relevant since we do not have the correct number 
of large space dimensions. Nevertheless, it is interesting to observe that a small coordinate can grow large 
even in the absence of any driving pressure. 



3.1.2 Ho^O, w = case 

Now, we investigate a universe filled with dust (w = 0) and an antisymmetric tensor potential {Hq ^ 0). 



Under this conditions, we substitute the ansatz (3.17) on equations (3.12), (3.13) to leading order as 



f> 



' ^^ ' U^ 



\V 



'0 -2dA 



A = (^A + H^e 



2„-2dX 



Using ansatz (3.17), (3.18), we obtain 



"P = 
A = 



2 + -jlogt + B 
1, 1, (Hl<f 



(3.23) 
(3.24) 

(3.25) 
(3.26) 



This analytic solution is plotted as gray colored straight lines in figure [3] for different values of Hq . Notice 
that Hq fixes the initital value of A(t = 1) in these solutions. 
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Figure 3: behavior of A and Lp when there is only the two-form field flux present {W 
Ho =0.001, Ha = 0.0001 ). 



0, K ^0, Hq^ 0, 



In this case we find that Hq is able, by itself, to induce the growth of a large scale factor, as can be seen 



in figure 3(a) In this figure we can see how the two- form field flux induces decompactification for different 
values of Hq. Since the two- form field flux is along 3 spacial dimensions, in the absence of winding and 
momentum modes, this field alone is able to induce the growth of 3 large spacial dimensions. We also notice 
that the moment in which the scale factor is able to " escape" the constant solution depends on the value of 
Hq. For larger values of it, the scale factor begins to increase earlier. 

This kind of scenario, in which the two-form field flux happens to be the dominant term, can occur if 
the pressure coming from the winding and momentum modes becomes negligible (P\ « 0) . This happens in 



generic situations, for example, when the scale factor remains near the self-dual radius, the dilaton goes to 
weak coupling or when the winding and momentum modes have annihilated. 

3.1.3 Ho^Q, w^Q case 

Finally we investigate the generic case when both the flux and the matter pressure are present. In order to 
find a solution when the antisymmetric tensor potential is present and the pressure fulfills the equation of 



state Px — wE, we use (3.121, (3.13) and the ansatz (3.141. Keeping only up to quadratic terms, we find 



A 
C 



2^2 

AC 

1^ 






ff2 

-"O ^-2d(Clogt+D) 



jj2^-2diC\ogt+D) ^ 



A^ 
1^ 



dC _ rT2-2d(C log t+D) 
t2 ^0^ 



By supposing that C = -5, we can eliminate the t ^ dependence on the equations. We obtain then 

^A 



1 
2d 



ff2 

-"0 -2dD 

2 



A2 

2 
1 _ A 
^d ^ ^ 
Substituting d — 5 explicitly and solving for A and D, we find 

5 — 3w 

A 



Hle-^'^^ 



'^{A^-\-Hle-^^^ 



H^e 



2„-6D 



3(w - 1) 
2 - 12w + 6u;2 
9(u;-l)2 



In this way we find a solution 



A 



5 — 3w 
3{w~l) 



-logt- - lo; 

3 



\ogt + B 

2- 12w + 6u;2 



We observe, on equation (3.32) that w is constrained by the inequality 

w <1- J-^ 0.1835 



(3.27) 
(3.28) 

(3.29) 
(3.30) 

(3.31) 
(3.32) 

(3.33) 
(3.34) 

(3.35) 



It is not consistent with w = 1/3 (radiation). This problem indicates that we cannot smoothly connect this 
solutions to those with Hq = 0. 



3.2 Perturbative solutions 
3.2.1 Ho^O, w^O case 

Due to the difficulty we just encountered, we construct perturbative solutions with non- vanishing flux starting 
from those with no flux. Using the solutions we have obtained for the case when Hq = 0, we treat the potential 
term due to Hq 7^ as a perturbation to the equations of motion. The small expansion parameter is 



m 



(3.36) 
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We expand the solution in terms of the small parameter S 



If ^ ipo + 5tpi + S'^tp2 + 
A = \o + 5Xi+ 5^X2 + ■ 



(3.37) 
(3.38) 



and substitute (3.371, (3.38) into the equations of motion. They describe perturbations around the solutions 



ipo, Ao obtained in (3.19), (3.20) 



From the power series expansion of the equation of motion (3.15 ), we get the differential equation for the 
first order terms in S 



<Pi = 'faV'i + dXo^i + -e 



-2rfAo 



W. 



Ai = ^oAi + Ao(^i + w{ipoipi - AoAi) + de -""^^ (1 - ^) 



After substituting ipo, Aq, <yJo, Ao into the equation, we obtain 

2 



^1 



1 + dw"^ 1 ^ '^'"'^ 



2dw 1 • 1 idiu ,„ 

l + dw^ ^2 



2 , • 2w 



1 + dw^ 



i-Vi- 



2dw 
1 + dw'^ 



t-'M) 



+rfriT^e-2'^^(l--) 
For Ai, we get a second order differential equation in this way 



Ai + 2t-iAi - driT^e-2'^^(l - ^) = 
We can integrate this equation easily. Defining a; = Ai, i = Ai we get 

i + 2t-^x = r^T^e-^'^^d{l - -) 
This is a differential equation of the form x(t) + f{t)x{t) — g{t) and the solution is given by 

J dtg{t)elf (')''' + c 

with a constant c. 

After the integration, we find two different class of solutions: 

• T^f^ 7^ 3 case. 

o^m , W^, Adw , 1 JdTii I 1 



^-"-'^"•'-IX'-T^)-'^ 



Adw 



^)-iri+<i™^+'-ct-i+co 



(3.39) 
(3.40) 

(3.41) 
(3.42) 
(3.43) 
(3.44) 

(3.45) 



(3.46) 



(3.47) 



1 + dw^ ' 

Here, the leading perturbation contains two different time dependent terms. For the perturbation to 
be small, the exponent on the first term should fulfill the condition 

\dw 



1 + dw2 



2<0 



(3. 
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If it is the case, the influence of the two-form flux induced potential is neghgible in comparison to the 
pressure of the string momentum modes. On the other hand, this condition is not satisfied for w = 
(pressureless dust) case where the perturbation grows as t^. In such a situation the solution Aq is 
unstable and the universe is decompactifled due to the presence of the two-form fleld flux. 

^'^^■2 — 3. This is the case for d — 3 and a universe filled with radiation {w = 1/d). 

(3.49) 







iogi c 


dX = 


s + c 


Iogi , 
s ^ +C0, 


t 



When t ^- cx), we find the leading perturbation as 5\ ^ 0{t^^) . Therefore the correction to the 
unperturbed solution is negligible at late time. 

3.2.2 Ho 7^ case with both momentum and winding modes 



As we observe in equation (3.10), the pressure coming from the winding modes and the momentum modes is 
multiplied by e"^. If |iJo| ^ 1 and \ip\ « 1, the scale factor experiences oscillations in the presence of winding 
and momentum modes. As (^ goes to weak coupling, oscillations stop and the pressure terms become small 
with respect to the Hq potential term. Before this terms becomes significant, the solution is characterized 
as 

93 «0, AsaO, A«0,e^< 1 (3.51) 

We define a small parameter 

^-f (3.52) 

Under this approximation, keeping terms to the lowest nontrivial order, we get 

(p = ^((yi^ + dA^ + e) (3.53) 

A = <^A + eA (3.54) 

We expand the solution in terms of the small parameter e 

(p = <y9o + e<^i + e^<^2 + • •• (3.55) 

A = \a + e\i + e^X2 + --- (3.56) 



After substituting (3.551, (3.561 into the equation of motion, we have differential equations at each order of 



the perturbation parameter 



V'o = 


1 -2 lj;2 
= i^vl + i^dXl 


Ao = 


= <^oAo 


^i = 


= (^o"^! + dAoAi -1- - 


Ai = 


= (^oAi -f Ao^Ji -f Ao 



(3.57) 
(3.58) 
(3.59) 
(3.60) 



12 



where the solutions for (po, Ao is given in equation (3.191, (3.201. Substituting them in (3.59), (3.60), we 
obtain 

2 1 

Ai - --Ai + co (3.62) 

with Co a constant. These equations are linear differential equations in ipi and Ai respectively. They can be 
solved by multiplying them by the integrating factor e'^J '^^t) — e(2iogt) _ ^2^ rpj^g solutions are 

Ai(t) = -^t + [const.) Xi{t) = -ft^ + {co7ist.)r^ + const. (3.63) 

Vi{t) = -t + (const.) ipi(t) — —-t^ + {const.)t^^ + const. (3.64) 

6 12 

We observe the following features: the solution ipo, Xq becomes unstable if we perturbed it with non- 
vanishing Hq. To leading order the solution in this regime behaves like ~ t^. This instability initiates an 
accelerated expansion of a universe away from the oscillating phase around the self-dual radius. However 
we also observe that the perturbation also affects the dilaton. As the perturbation becomes dominant, the 
dilaton begins to grow and goes to strong coupling. This indicates that a bounce on the dilaton has been 
produced. A result like this looks problematic, since a bounce on the dilaton leads to a violation on the 
positive energy condition as was noted in |16| . We may not be able to trust our solution there as it also 
takes the dilaton to strong coupling. This behavior can be observed directly in a numerical solution of the 
equations (figure HI). We begin with a string gas of equal number of winding modes and momentum modes. 
Before the dilaton goes to weak coupling, the scale factor oscillates around the self-dual radius. Once the 
dilaton reaches weak coupling region, the oscillations stop and the scale factor stabilizes. Then the potential 
induced by the two-form field flux becomes dominant and the solution begins to grow as predicted by the 
perturbed solution. In the next sub-section, we investigate the effects of the string interaction on these 
problems through the Boltzmann equations. 

3.3 Effects of string interactions 

Up to this moment, we have considered situations in which the winding and momentum numbers are frozen 
at their initial values. When the string gas falls out of equilibrium in an expanding universe, winding strings 
in the gas can interact and begin to annihilate. In this section we incorporate, together with the two-form 
field flux induced potential, the Boltzmann equations that take account of the interaction among strings. 
These equations, derived by Polchinski [27], are shown below 

W = {W^-{W)^) (3.65) 

K = (if2_(i^)2) (366) 

TT 



We combine these equations with (3.9), (3.10) and evolve the system numerically. The universe we consider 



is filled with gas of strings that begins at the self dual radius with equal initial winding and momentum 
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Figure 4: behavior of A in presence of U{X), winding and momentum modes {W ^ 0, K ^ 0, Hq 
Ho = 0.01, Ho = 0.001 ). 



numbers {K = W). The initial conditions are A ~ 1, (/i ~ — 1 and the dilaton is going from strong coupling to 
weak coupling. The numerical results including the effects of the Boltzmann equations are presented in figure 
Isj We have plotted the behavior of the scale factor \{t), the dilaton (p{t), the winding number w{t) and the 
momentum number k{t). In figureplwe observe that, as A(i) grows, the winding modes begin to annihilate. 
Then, there is not enough pressure to make the universe contract and experience bounces. Instead, the 



contribution from the two-form field becomes dominant and the scale factor tends to the solution (3.251 



(3.26) with vanishing pressure where the scale factor grows large due to the flux induced potential. 



The behavior of the winding and the momentum number is as expected from the following characteristics 



of the Boltzmann equations (3.65), (3.66). As the dilaton goes to weak coupling, the interaction rate goes to 



zero and the values of the winding and momentum numbers become constant. When the scale factor grows 
large, winding modes annihilate more efficiently because their interaction rate goes as the exponential of the 
scale factor. On the contrary, the rate of annihilation of the momentum modes becomes smaller because the 
interaction rate between them decays exponentially with the scale factor. In fact this asymmetry between 



winding and momentum modes can be observed in figure 5(c 



The result obtained by taking account of the Boltzmann equations suggests an interesting scenario when 
homogeneous H^^\ is present. If the winding modes annihilate rapidly enough, the effect of 2-from field 
flux becomes important even at early times. The annihilation of the winding modes could take place even 
in a loitering phase. In that case the two-form field flux becomes dominant and the the expansion of three 
large spacial dimensions is realized. We emphasize that this mechanism is different from Brandenberger-Vafa 
mechanism as the presence of homogeneous 2-form field flux is crucial for three spacial dimensions to grow. 



Without it, the universe remains to be of microscopic size as the blue line in 5(a) indicates 
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Figure 5: (a), (b) behavior of A and ip with initial condition K = W ioi Hq = Q, Hq = 0.001 and Hq — 0.1. 
The effect of the Boltzmann equations is included, (c) Evolution of the winding number and momentum 
number is plotted (the solutions overlap for the three cases considered). 
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4 Effects of constant B 



flV 



So far, the effect of the two-form field has entered only as a modification to the usual dilaton gravity action, 
as in [Tl]. The string gas model, as it stands, couples the modified action of dilaton gravity with that of a 
gas of strings. In this approach the effect of the background field B^^, over the string spectrum is usually 
neglected. The correction on the energy of the string goes as 0{B), thus, this approximation is valid for 
weak fields. In dilaton-gravity, the contribution of Bfj,^ to the action enters via U{X) ~ \H\'^ = |di?p. In this 
case, even if B remains small, H is not necessary so as the space-time variation of B^^^, could be large. 

In principle, if we know the two-form field in terms of the scale factors, we can determine \H\'^ as well 
as their effect on the string spectrum. Wc can then make use of the adiabatic approximation to study the 
time dependence of the compactification radii and get the equations of motion. In practice, a homogeneous 
solution for supergravity is given in terms of H. This presents a problem since we need Bfj^,^, not H, in order 
to get the string spectrum. 

With this prospect, we investigate the simplest case, that of a constant i?^^. In this case, the H dependent 
term on the supergravity action vanishes as well as the contribution to the equations of motion. Nevertheless, 
since strings carry charge under the gauge field, the effect of B field on closed strings wrapping the compact 
dimensions is non-trivial. 

The Polyakov action in the presence of an antisymmetric field B^i, 



S = - 



1 



d^aldaX'^d'^X^ ~ e''''B^,daX^dbXf'] 



(4.1) 



yields the equations of motion 



{dl - dl)X^\cj,T) = --H'^^.e-'daX^dbX-' 



(4.2) 



with Hf,^x = df_,Bi,p 
equation 



dpBni, 



dijBp^ . Then, for a constant B^u we obtain the usual two dimensional wave 



{dl - dl)X'' ^ 
that allows us to give the solution in a Fourier-Laurent expansion 



X'^' = 



— {aa + ao) + 



(ao - ao) + 






e Mr+'y) _^ ^l^^~in{T-a)^ 



i(r+cr) 



^Mg-™(r-<T)N 



(4.3) 

(4.4) 
(4.5) 



^(a„e 

As it turns out the zero-modes are the only ones that are affected by the B^^ field. The components of the 
energy momentum tensor and their zero modes are given by 

1 



rpQl _ 



2na' 



1 



a'pt" 



B^'jiwRy 



tE«^ 



l(T+cr) 



-Hat^e" 



^{T—a■)^ 



(wR), 



n/0 



y^nfi^ 



n/0 



— in(r+cr) __ 



-iniT — a) 



'^nfi^ 



(4.6) 
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V-^ jzcro modes — ^ '^i \ ^- / jX^n ' ^—n ^ ^ — n ' ^n ^ — n ' ^n ^n ' ^—n) 



(4.7) 



T — 



1 



1 



-{X -X + X' ■ X') 



a'pf'-Bf'jiwRy 






-»n(r+a)^ Ug-^„(r-^)- 



(u;i?)'' + V T E«^""^'^'^ - <e-™(^-'^): 



n/0 



2n 



(4.8) 



V-*^ yzcio mode 



1 



- (ay - B°j{wRyf + (aV - B'j{wRyf + (wi?)*(wi?,) 






(4.9) 



Imposing the physical constraint that the energy momentum tensor must vanish, we get the energy spectrum 
for the string 



_ B°wm 1 



and the level matching condition 



N -N ^ n'w, 



(4.11) 



Since all the spacial coordinates are compactified with radius R^ , the momentum is quantized as p* = (n/RY , 
where i denotes the spacial index. 

In order to be able to solve the equations of motion, we need to assume some initial winding and 
momentum distribution of the string gas. The constant B^^, field could be either electric or magnetic type. 
We find that the effect of electric type field is very interesting as there is a critical value for which the string 
tension vanishes for winding modes. 

4.1 Constant electric type field 

Let's consider the case of a homogeneous electric type field in 3-spacial dimensions, with B = Bqi = i?02 = 
_Bo3- In order to demonstrate the most dramatic effect, we assume that strings are aligned in the direction of 



the electric type field. If this is the case, from (4.10), the energy for the winding modes in (2.21 ) is modified 
as 



Ew = 2d(l - V3B)We^ 



(4.12) 



with W the winding number and d ~ 3. From this equation, we see immediately that the effect of the field 
S^i, is to reduce the energy of the winding strings. Also, it follows that B is constrained to take values 

1 



0< B < 



(4.13) 
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Figure 6: behavior of the dilaton and scale factor for 5 = 0, S = 0.5, B — I/a/S (critical electric field) and 
W ^K. 

In particular, when the inequality is saturated B — l/-\/3, the energy of the winding modes vanishes. As 
the pressure they exert also vanishes, the spacial dimensions are expected to expand freely because of the 
presence of the momentum modes. 

In figure |6] we have plotted the numerical solution for different values of B without including the effect of 



the Boltzmann equations (3.65). For vanishing B the momentum and winding modes make the scale factor 



oscillate around the self-dual radius. With small S ^ 0, the solutions oscillate around positive values of A. 
As we get closer to the critical i3, the solution bounces and then stabilizes. When we reach the critical value 
B — l/vo, the pressure from the winding modes becomes zero and A expands just like a universe filled with 
radiation (momentum modes). 

5 Conclusions 

In this work we have investigated some effects of the introduction of a two-form field into the model proposed 
in [TB] . This model provides a bouncing and cycling cosmology and also the possibility of long loitering phases. 
It avoids singularities at finite times but fails to realize three large spacial dimensions from Brandcnbcrger- 
Vafa mechanism. Having this in mind, we have included a two-form field into the action, since it may provide 
an alternative mechanism for the decompactification of 3 spacial dimensions. We have considered two cases: 
homogeneous flux B^^yp and constant gauge field -B^j/. 

5.1 Homogeneous Hp^^p 

In order to make the model compatible with T-duality, as the string gas model requires, we have adopted 
a phenomenological modification on the potential induced by the two-form field entering the gravity action. 
The modified potential is non-singular at A = and reduces to the correct one when the scale factor |A| is 
large. In addition, it provides a repulsive potential that can make the universe expand. 
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In the investigation of the behavior of the scale factor and the dilaton under the influence of the two-form 
field flux, we find two different cases: 

• Matter dominance: 

At early times the scale factor can experience bounces as it is governed by the presence of winding 
and momentum modes. In section 13.2.21 we have observed that the effect of the two form field is not 
significant at this early stage of the universe. If we assume only the presence of the momentum modes, 
the late time solutions reduce to those already found in dilaton cosmology. If this solution is perturbed 
by the introduction of the two-form field flux potential, its influence vanishes as t — >■ oo. On that 
account, this kind of solution is stable under the perturbation and the effect of Hf^,jp is negligible as 
the universe expands. 

• Two-form potential dominance: 

We have obtained the late time analytic solution for vanishing matter pressure and non- vanishing H^^p. 
This solution corresponds to an expanding universe, where the initiation time of the expansion is set 
by the parameter Hq. This analytical solution matches the leading behavior of the numerical solution 
for the equations of motion. 

In generic situations the contribution of the matter pressure becomes negligible and the scale factor 
becomes constant. This occurs when the dilaton goes to weak coupling, the oscillations on the scale 
factor stop or the expansion of the universe comes to a halt. Such a possibility is enhanced if we 
consider the effect of interactions between strings. As momentum and winding modes can annihilate, 
it drives the pressure to vanish. In all of the above cases, the effect of the matter pressure vanishes and 
the scale factor becomes approximately constant. Introducing a Hf^i^i^p flux, we find that the constant 
scale factor solution eventually becomes unstable and the scale factor begins to grow as 

\^e (5.1) 

This kind of scenario occurs whenever the dilaton goes to weak coupling and the scale factor settles to 
a constant value. This behavior is remarkable, since it produces an accelerated expansion analogous 
to the inflationary universe. However we also need to address the issue that the perturbation to the 
dilaton also goes as t^ . Thus the dilaton may eventually bounce and go to strong coupling. The string 
interaction effect through the Boltzmann equation is observed to resolve these problems as in figure [5] 

5.2 Constant B^^ 

We have also investigated the case of a constant Bp^, in order to test how its presence affects the action of 
the string gas. For a constant field, the equations of motion for srings reduce to the usual one without Bp^. 
It is straightforward to include the effects of a constant B field by calculating the spectrum of the string. 
The inclusion of a constant Bp^, has some interesting consequences, one of these is that there is a critical 
value which makes the energy of the winding modes aligned with Bp^, field vanish. 

We had expected the modification induced by Bpi, to be significant since its presence makes the energy 
and the pressure of the winding modes vanish at a critical value. In fact our numerical results indicate that 
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the behavior of the scale factor could be significantly affected. The spacial directions expand like radiation 
dominated universe even with the presence of the both momentum and winding modes. 
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A Relation between string frame and Einstein frame 

In this appendix, we summarize the relation between string frame and Einstein frame in our setup. The 
parameter d in this appendix which counts the number of the spacial dimensions should be put rf = 9 in 



superstring. From the conformal transformation (2.7) and the corresponding metrics the relation between 
scale factors is 

a = -^ + A, /3 = -/^ + ^ (A.1) 

dP^e^^dt^ (A.2) 

Also, the shifted dilaton is defined by 

d 

(p = 0-^A, (A.3) 

1=1 
Because of the presence of Hfj^^x for the superstring in ten dimensions, the spacial coordinates factorize as 
T^xT^. Defining A, = A for i = {1, 2, 3} and A^ = zy for j = {4, • • • ,9}, we have if = 4)-3X-6iy. Using the 
Einstein equations and the solution for the homogeneous two-form field, we obtain the equations of motion 
for the superstring case {d — 9) 



^^~3X^-6u^ = 


e'^E + UoiX) 








(A.4) 


if-3X^- 6i>^ = 


\e^E 








(A.5) 


A — ipX = 


ipA + f/o(A) 








(A.6) 


i) — ipi/ — 


2' ^'^ 








(A.7) 


with 












C/„(A) ^ ^Hl,^ = ^H'^e-'^ 








(A.8) 


Accordingly, the equations of motion in the Einstein fr 


ame are 










fin^'2 1 qfin''^' 1 qn'^'2 — "^ 1 r''''^--!^ 


3"-6^iJ+ii/2g- 


6-3^ 






(A.9) 


a — 1 


a" + a' {3a' + 6(3') = +|i/2e 6" ^"-i 

8 


+ 1 p{^^-3a 

2-8' 


"^^(-5£; + 5Pa- 


-6P,) 


(A.IO) 


/3" + /3'(3a' + 6/3') = -^i/^e ^^ .% 

8 


-+- p(d-l)f 3a 

2-8' 


"^''(S^;- 


3Pa + 2P,) 


(A.ll) 


0" + </.'(3a' + 6/3') = H^,e «" - ^ + - 


^e<-')*-^"-«^(3P,-f6P.) 






(A.12) 



20 



In the string frame, the equation of motion for A contains the dilaton and its time derivative but it does not 
contain v terms. In the same way the equation of motion for v is independent of A or its time derivatives. 
Then, the equations of motion for the scale factors v and A decouple and we can proceed to solve them 
numerically. In comparison, in the Einstein frame, the presence of li^lv\ makes the scale factors couple to 
each other. In spite of this unfavourable characteristic, the equations of motion in the Einstein frame are 
also useful, both when trying to solve the equations of motion and also for clarifying the interpretation of 
the solutions. 

In the Einstein frame the field U^^^x is included in the equation of motion for both a and /3. By looking 



at the sign of the Ho term in (A. 10 1 and (A. 11 1 we can see that the two-form field induces an anisotropic 
expansion on the scale factors, with and a being driven towards positive values while /3 goes towards 
negative values. Also, while in the string frame it is possible to find solutions to the equations of motion in 
which V becomes constant, this does not imply that the physical scale factor is fixed because it remains to 
stabilize the value of the dilaton. This can be seen directly from the relations of the Einstein frame to the 
string frame, where, in the case oi v — constant we have 

/3 = — -^ + const. (A.13) 

a— 1 

That is, unless both the dilaton </> and v are constant in the string frame, there is no solution with /? = 
constant in the Einstein frame. 
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